EXPECTED VOLUME OF INTERSECTION OF WIENER 
SAUSAGES AND HEAT KERNEL NORMS ON COMPACT 
RIEMANNIAN MANIFOLDS WITH BOUNDARY 



M. VAN DEN BERG AND P. GILKEY 

Abstract. Estimates are obtained for the expected volume of intersection of 
independent Wiener sausages in Euclidean space in the small time limit. The 
asymptotic behaviour of the weighted diagonal heat kernel norm on compact 
Riemannian manifolds with smooth boundary is obtained in the small time 
limit. 



1. Introduction 

Let Bi, Bk be k independent Brownian bridges in R m associated with the 
parabolic operator — A + J^, and with -Bi(O) = Bi(t) = for 1 < i < k. Let if be a 
compact, non-polar set in R m . For 1 < i < k, let 

Wk(t) = {Bi(s) + y.0<s<t,yeK} 

be the corresponding pinned Wiener sausages. These random sets are Borel mea- 
surable with probability one, and we let 

Z k , m (t) := E 1 ® ... ® E k 

It is well known that Zk, m (t) is related to the virial coefficients of a quantum 
system of obstacles K at inverse temperature t. For example, G. E. Uhlenbeck [T5] 
calculated the asymptotic behaviour of Zk im {t) as t —* oo in the special case where 
k = 1, m = 3, and if is a ball. I. McGillivray ^1] obtained the full asymptotic 
series as t — > oo for k — 1, m > 3, and K an arbitrary, non-polar compact set. The 
special case where k = 1, m = 2, and if a ball was studied for t — * oo by M. van 
den Berg and E. Bolthausen 0]. 

In this paper, we analyse the asymptotic behaviour of Zk yTn (t) as t — > 0. Let 
PK m (', S ') denote the heat kernel for R m given by 

p Rm (x,y;i) = (4^r m / 2 e-l^l 2 /( 4t > , 

and let pr™_ k denote the Dirichlet heat kernel for the open set R m — K. By the 
Feynman-Kac formula, we have that for x € M. m — K, 

(La) [B(s) nif = 0,O<s<<] = {47rt) m / 2 p Rm _ K (x, x; t) , 

where B(s), < s < t is a Brownian bridge with £?(0) = B(t) = x. It is easily seen 
that 

(Lb) < p K m._ K (x,x;t) < p Rm (x, x; t) , 

and that 

(l.c) Z k , m (t) - {4nt) km / 2 [ {p Rm (x,x- 1 t)-p Rm ^ K {x,x;t)} k dx , 

where the Dirichlet heat kernel is extended to all of R m by putting 
Pr™~k{x, y; t) := for x 6 K and or y <E K . 
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This extends the equality of formula (|l.a|l to all of R m . 

It follows from the results in jS] that if k = 1, m > 2, and K is compact with 
smooth boundary, then there exists an asymptotic series as t — > 0: 

3=0 

where the coefficients cij have been computed for j = 0, 4. For example, 

ci.o = / Idx, ci,i = [ Idy, and ci i2 = \ [ L aa dy , 
J k 1 J ok o Ja/f 

where (iy is the surface measure on dK , L aa is the trace of the second fundamental 

form on dK, and where dK is oriented with an outward orientation, i.e. the inward 

normal is chosen on d(M — K). 

Let f (t) be defined for t > 0. We say that £(*) - £„> £„i™/ 2 as i ^ if 

.7 

f (t) = £ £„/;"/ 2 + 0(< (J+1)/2 ) for any J G N . 

n=0 

The main results of this paper are the following two theorems: 

Theorem 1.1. Let K be a compact set in R m (m > 2) with smooth boundary 
dK. Let k G N be arbitrary. Then one has that Zf. <m (t) < oo for all t > 0, and 

c k , =j K ldx, CkA = \(l) 1/2 j^ldy, 
ci,2 =~\ ( L aa dy, 
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(5/4 



anrf /or A; > 2 

Cfe,2 



1 



2fc 2 V4(fc- l) 3 / 2 

1 arctan(/c — 1)~ 1//2 



L aa dy. 

dD 



2k(k-l) 2{k- If/ 2 

Theorem 1.2. Let M be a compact m- dimensional Riemannian manifold with 
smooth boundary dM . LetpAi(-, ■;■) denote the heat kernel for the Laplace- Beltrami 
operator acting on L 2 (M) with Dirichlet boundary conditions on dM. Let f be 
smooth on M and let k G N. 

(1) Fort^O, / {p M {x.x-t)) k f{x)dx ^ (^t)- km ' 2 Y^a k ^/ 2 . 

(2) a k ,a = / fdx. 



M 

k 



(3) «M = ^£(-l)^ 1/2 ( * ) / 

1=1 \ ' ^ 9 

k f 

(4) a fc! 2 = - / frdx 

6 Ja/ 

-\t^i f (1)d y + {- k e + \t 

2 1=1 J aM I o ^ 



/di/- 

dM 



(-1) 



(fc ^(4^/2 2^+1) C 2f3/2 ) + JIL/ 71 " 1 ^ 
where r is £/ie scalar curvature, and where /W i/ie normal derivative of f with 
respect to the inward unit normal vectorfield on dM. 
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It is possible to obtain the formulae for the Ck j,j > 1 in Theorem 11.11 directly 
from the formulae for the ak,ji3 > 1 m Theorem 11.21 by noting that the formulae 
in Theorem 11.21 also hold for non-compact m-dimensional Riemannian manifolds 
provided DM is compact and smooth and / is smooth and has compact support. 
Specializing to the case where M = M™ — K, f compactly supported and identically 
equal to 1 in a neighbourhood of dK and noting that pm_™ (x, x; t) = (47rf) -m / 2 yields 
a closed set of equations relating the Ckj and a^j by equating powers of t: 



(l.d) 



3 e 



This paper is organized as follows. In Section we consider a half space and a 
suitable localizing function /. The formulae for a,k,o and au,i in Theorem 11.21 follow 
from Lemma 12. II One of the two boundary terms in a,k,2 also follows from Lemma 

i2~n 

In Section we study a planar region in Euclidean space. In Lemma 13.31 we 
determine the coefficient of J dD L aa dy in Cfc 2 - This completes the proof of Theorem 
11.11 The second boundary term in 0^,2 is determined in Lemma 13.41 by computing 
the coefficient of J gM fL aa dy. This completes the computation of the boundary 
term in 0^,2- It would be a tedious combinatorial exercise to show that these 
coefficients are consistent with formula <|l.d[) . 

In Section 0] we complete the proof of Theorem 11.21 by determining the interior 
term in au,2- This uses the results of McKean and Singer JSj- The arguments used 
in this section together with existing formulae [7] would suffice to determine dk,e 
for all k and for I < 6 if the boundary of M is empty. The crucial difficulty in 
this paper, however, is the determination of the boundary terms. These are not 
accessible by the methods used in Section 0] 

In Section El we obtain results for the expected volume of intersection of Wiener 
sausages which are not pinned, and which complement Theorem ll.il 



2. Half space calculations 
Let Pk+ ( x ,x;t) be the Dirichlet heat kernel for the positive half line R + . 

Lemma 2.1. Let k G N, and let f be smooth. For t — > 0, we have 
/ {PR+(x,x;t)} k f(x)dx 
- (47rf)- fe / 2 f:.r:d.r ■ -/ < 1 | 



where 



1 / -,\P f k 

2 ^ 

1=1 



{-if ( \ ) r (0 r(j)-H-^ 
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Proof. We have that 

p{x,x;t) = {A-Kt)- x / i {l-e- x *l% 
/(*)~Sy W) (0y as x^O, 

3=0 3 ' 

(1 - er'^l'f = 1 + Y,(-lf ( * ) e- fe2 /*, 

Lemma f2 . II follows from these identities. □ 
Note that by identity (1.45) of [HI 

3. Computations for planar regions and proof of Theorem 11.11 

Lemma 3.1. Let K be a compact set in R m , and let Zk, m (t) be given by Ijl.cjl . 
Then Zk, m (t) < oo for all t > 0. 

Proof. Since p^m (x,y : t) < (Ant)-™/ 2 we have by (|T7b)l . (fTcj) that 

Zk. m (t) < (AiTt) m/2 / {p Rm (a;,a:;t) -p Rm _K(a;,^;0} rfa; 

JR™ 

jR m Jr™ 

^ 3 - a ) < 2 (2+m)/2 f I {mm ^y. t /2)-p Rm _ K ( X ,y;t/2)}dxdy 

Jl™ JR m 

= 2< 2+m )/ 2 / (l- / pRm-jfCa, y;t/2)d!/}d!B 

JM m -K I JWL™-K J 
+ 2 (2+™ ) /2| K | j 

where we have used the semigroup property of the heat kernels for R m and M m — K 
respectively. The first term in the right hand side of l|3.a[l is the amount of heat in 
R m - K at time t/2 if K is kept at fixed temperature 2( 2+m )/ 2 , and if K m - K has 
initial temperature 0. This term is finite by the results of F. Spitzer [T2|. □ 

A key ingredient in the proofs of Theorems 11.11 and 11.21 is an estimate for the 
Dirichlet heat kernel obtained by R. Lang [§] and H. R. Lerche and D. Siegmund 
|llj . For related results see also ^21- Before we state their results (Theorem 13.2(1 
we introduce some further notation. For an open, bounded and connected set D in 
R 2 we define the distance function S : D — > (0, oo) by 8(x) = min{|x— z\ : z £ dD}. 
Suppose dD is of class C 3 . Then there exists Ed > such that for all x £ D with 
5{x) < Ed there is a unique s £ dD, depending on x, with \x — s\ — 5{x). Hence 
for all such x £ D there is a smooth parametrization x — > (s, 5) where s £ dD is 
parametrized by arc length. 

Theorem 3.2. Let D be an open, bounded and connected set in M 2 with C 3 bound- 
ary. Let x £ D be such that 5(x) < Ed/2. Then as t — > 
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p D (x, x; t) = (Airt)- 1 { 1 - e- s ^ ' l - L aa (s)S(x) 2 r^ 2 / e"" dr, 

[ Js{x)t~ 1 / 2 J 

(3.b) +0(1), 

where the remainder 0(1) is uniform on {x 6 D : d(x) < Ed/2}, L aa is the curva- 
ture of dD at s, and where dD is oriented by an inward unit vector field. 

Lemma 3.3. Let D be as in Theorem iS.'A and let k £ N. Then as t — ► 

2 

(3.c) (Ant) k / {p M 2(i, i; «) - i; t)}^! = V c fcj ( ,/2 + o(i), 

where c k ,o = \D\,c k ,i = |(f ) 1/2 J azJ Ids, 

ci,2 = -g / L aa (s)ds, 

JdD 

and for k > 2 

_ J 1 fc / 7T 

° M t~ 2fc + 2 v 4(fc-l)3/2 

1 arctan(fc — l) -1 / 2 



2fc(fc-l) 2(fc- 1)3/2 J j J 9D Laa{s)ds - 

Proof. Let e e (0, |) and denote D £ = {i e D : i5(i) < i £ }. Since 

(3.d) Pm2{x,x;t) > PD {x,x;t)>p R 2(x,x;t) - {nt)- 1 e~ 5 ^ 2 ^ 

(see Theorem 1 and estimate (ii) in pQ) we that for e < \ and t — > 

(3.e) / {p R 2(a;,a:;t) -p D (a;,a;;0} fe c?a; = 0(e^* 2 ^ 1/3 ). 

Jd\D s 

It remains to calculate the contribution from the set D e to the integral in the left 
hand side of l|3.cf) . Since the 0(1) remainder estimate is uniform in x we have that 
this remainder estimate contributes at most 0(t 1+£ ) to the integral in 1)3. cfl . Hence 
for t sufficiently small we have by (|3.e|l 



Z k , 2 (t)= I da [ drie-r 2 ' 1 

JdD JO I 



+ L aa (s)rH~ 1/2 / e-^dr, (1 - rL aa {s)) + 0(< 1+e ). 

Jrt- 1 / 2 J 

Since /„* drr 4 ^ 1 = 0(t 5E_1 ) we have that for e 6 (fil) the contribution from 
such a term is o(i). Hence for t — » 

rt 6 r /-oo 



rt -l/2 



Z*, 2 (*) = f ds f dAe-^l' + ke-^-^/'LMrH- 1 " 

JdD Jo I 

(3.f) x(l-ri aa ( S )) + 0(i 5e - 1 ). 

Furthermore for e < 1/2 and t — ► 

- / ds / dre~ kr2/t rL aa {s) ~ -— / dsL aa {s), 

JdD Jo iK JdD 

dsL aa (s) f drrH- 1 ' 2 f e"^d»j ~ | / dsi aa (s), 

OD JO Jrt- 1 / 2 JOD 
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and for k > 2 



2 



dsL aa (s) / drfee-^-^/Vr 1 / 2 / 'd V 

dD JO Jrt- 1 / 2 

§ / dsL aa (s) H drke-W^t- 1 '* r e-^drj 

dD JO Jrt- 1 / 2 



dsL aa (s) I drtiij 2 + k - iy 2 

dD Jl 



kt ( 7r 1 arctan(fc — 1) 1 / 2 



2 \4(fc- 1)3/ 2 2k(fc-l) 2(fc-l) 3 / 2 j J 9D dsLaa{s) - 

It is easily seen that the term with L aa {s) 2 in (|3 - ff> contributes 0(i 3 / 2 ). Collecting 
the contributions from j3.g| - I3.hjl we arrive at the conclusion of Lemma 13.31 □ 

Proof of Theorem M.lX Since the contribution to the asymptotic series for Zk^ m (t) 
comes from & t e ,s < 1/2 neighbourhood of dK we can identify the coefficients 
Cfc.OjCfc.i f° r m — 2 by Lemma |3.3I By standard results of invariance theory they 
also hold for m > 2 [?]■ Moreover Zk,m(t) has a full asymptotic series. The proof 
of Theorem 1 1.1 1 is complete by Lemma f3. II □ 

Lemma 3.4. The coefficient of the term J dD fL aa dy in the expression for ak,2 in 
Theorem \l.£\ is given by 

k 1 >^ / k \ . p\ f n 1 arctan£ _1 / 2 \ 1 

Proof. It suffices to consider a planar region with / = 1 as in the proof of Lemma 
E2 Let e € (2/5, 1/2). Then by {3~Hj) 

dx{p D {x,x;t)} k ~ (Ant)- km / 2 \D - D B \. 

I D — D e 

The O(l) remainder in (|3.b(l contributes at most 0(t~ k+1+£ ) to the integral over 
the set D e . Hence for 2/5 < e < 1/2 

dx{p D (x,x;t)} k = {4irt)- k [ ds [ dr(l - rL aa {s)) 

D E JdD Jo 

x 1 - e" 1 " 2 /* _ r H-^L aa { S ) / e^dr, + 0(t- k+1+£ ) 

I Jrt- 1 / 2 J 

= (4^)- fc / ds f dr(l - ri aa (s))((l - e^*) k 

JdD JO I 

r°° "i 

Jrt- 1 / 2 J 

by estimates similar to the ones in the proof of Lemma l3~3l For 2/5 < e < 1/2 and 
t — > we have 

(4^)~ fc /" da / <fr(l - rL 00 (s)) = (47rf)- fc |£> £ |, 

JdD JO 



lnt)~ k I ds [ dr{(l 

JdD Jo 



(4irt)- k I ds I dr{(\ - e~ r '*)* - 1} 
2- 1 (4 7 rr fc (^) 1/2 E( ) )(-l)^ 1/2 , 
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and 

(3.i) -(47rf)- fc / dsL aa (s) I drr{(l - e^' t ) k - 1} 



dD 

k 



Furthermore 



2- 1 (47rt)- k tJ2 ( \ ) (-ly-'r 1 J^dsL aa (s). 

-A;(47rt)- fc / d Sj L aQ (s) / drr 2 ^ 1 / 2 / e~ r ' 2 
Jan Jo Jrt- 1 / 2 



ch) 



~ -^i(47rf)- fe i / dsL aa (s), 
6 Jan 

and so it remains to compute ( for k > 2) 

(3.j) -fc(47rf)- fc / dsi aa (s) / dr 

Jan Jo 



x 



-fc(4 7 ri)-' ! / dsL aa (s) / dr 



x 

fe-i 



r fc-i 

JdD l=l 



x(-iy 



1 arctani? 1 / 2 



,4£ 3 / 2 2£(*+l) 2£ 3 / 2 

Collecting all of the above we have shown that the formulae for a^o, Ofc,i and 
in Theorem 11.21 hold for planar regions with C 3 boundary. In particular we have 
shown that, by collecting the terms (|3-H3.j| ), the coefficient of the term J gD fL aa dy 
in the expression for ak,2 is given by Lemma 13.41 The scalar curvature term in 
the expression for au,2 in Theorem II .21 (A) will follow from the results in Section^ 
below. □ 



4. The interior terms 

In order to obtain the contributions from the interior of M to the coefficients 
Ofej- in Theorem ll.2l it suffices to consider Riemannian manifolds without boundary. 
The general setting is as follows. Let D be an operator of Laplace type on a smooth 
vector bundle V over a Riemannian manifold M without boundary. One can express 
D in terms of geometrical data in an invariant fashion as follows. There is a unique 
connection V on V and there is a unique endomorphism E of V so that one can 
express D in the Bochner formalism by writing: 

D = {;/"V,V, + E} . 

Let Pm{x, x\ i) be the heat kernel, and let F be a smooth endomorphism of V. The- 
orem generalizes to this setting to yield the existence of a complete asymptotic 
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expansion 

~ oo 

/ T± Vx {(p M (x,x;t)) k F(x)}dx~(4nt)- km / 2 J2a k!2j (F,D)ti, 

where the coefficients a k ,2j are locally computable. There are no half integer powers 
of t since the boundary of M is empty. Consequently a k ,2j+i = 0. 
If k = 1, then we have a complete asymptotic expansion 

oo 

(4.a) p M (x, x; t) ~ (47rf)" m / 2 £ e 2v {D)t v 

u=0 

where the e 2v are locally computable endomorphisms. Let r := Rijji be the scalar 
curvature. Then one has the following well known result (see, for example, the 
discussion in [7]): 

Lemma 4.1. Let D be a operator of Laplace type on a closed Riemannian manifold. 

(1) a hn (F,D) = J M Tr{Fe n (x,D))dx. 

(2) e Q (x,D)=id. 

(3) e 2 (x,D) = ±{rid+6£}. 

Remark 4.2. We note that similar formulae are available for e^ and ee [7|. 

Raising the asymptotic expansion in formula (|4.a() to the k th power and applying 
Lemma f4. II then yields immediately the following: 

Corollary 4.3. Let M be a closed Riemannian manifold. 

(1) a k , (F,D) =a lfi {F,D). 

(2) a k , 2 {F,D) = k ai , 2 (F,D). 



(3) a M (F, D) = ka 1A (F, D) + J Tt{F^^-{t 2 + \2tE + 36E 2 )}d. 



x. 



5. Expected volume of intersection of independent Wiener sausages 

In this section we obtain results analogous to Theorem 11.11 for the unpinned 
Wiener sausage. Let (3 1 (s) , . . . f3 k (s) , s > be independent Brownian motions in 
R m , and let K be a compact set in l m . The Wiener sausages S]^ (t), . . . , S^(t) are 
the random sets defined by 

S i K (t)={(3 i (s)+y: 0< S <(,|/e4 i = l,...,*. 

These random sets are compact and Borel measurable with probability 1. Define 

" k 

Q k ,m{t)=K®---®K 

It is well known that 



k 

Qk,m(t) = I dx\l- / p R m_ K (x,y : t)dy 

and that Q k)Tn {t) < Qi, m {t) < oo (CHI)- We note that Q k , m (t) has the following 
analytic interpretation. Let u : R m — K X (0, oo) — * R be the unique weak solution 
of Au = ^ with initial condition u(x; 0) = 0, x 6 R™ 1 — K and boundary condition 
u(x; t) = 1, x S dK , t > 0. We extend u to all of R m x (0, oo), and note that 



u(x;t) = l- pR m ^ K (x,y;t)dy. 

JR m -K 

Then 

Qk,m(t) = \\u(-,t)\\ k k . 

The large t behaviour of Qk,m(t) has been investigated for m = 2 in ^H] and for 
to > 3 in 0. 
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Theorem 5.1. Let K be a compact set in M m with C°° boundary. There exists an 
asymptotic series as t — > 

oo 

Qfe,m(t)-^&fe J t j/2 , 
3=0 

where 

(5. a) b kfi = / Idx, 

J K 

(5.b) 

6 fc . 1 = 2 2 +V 2 - fe )/ 2 r((fc + l)/2) / <%.../ ^ fe (r / 2 + ... + r ? 2 )-( fe + 1 )/ 2 / ldy, 
6 fc ,2 = -2 2+fe 7r (2 - fc)/2 (fc - 2)r((fc + 2)/2) 

(5.c) 



x 



dm... dmivl + ■■■ + vtr {k+2)/2 / L aa d y . 

1 Jl JdK 



Proof. It suffices to prove (|5.aH5.c|) for a disk in R 3 . The general case for compact 
K with C°° boundary follows by invariance theory. For the disk D 3 centred at the 
origin it is well known [£] that 

/>oo 

u(x\t) = 2ir- 1/2 r~ 1 / e^dr], r = \x\ > 1. 

J(r-l)/(2t 1 / 2 ) 

Hence for K = D 3 we have that 



(r-l)/(2t 1 /2) 



The main contribution to Qk.3(t),t — > comes from the interval [1, 1 + i e ], where 
< e < 1/2. Expanding r 2 ~ k near 1 yields 

Air / 9 \ ^ Z" 00 f/" 00 1 



™ /• OO /> CO 



2 



8tt 



>7r l/2 

0(i 3/2 ), 



/» OO /"CO />OC 

(fc-2) / drr k+1 d m ... dn k e- r2 ^l+-+ti) t 
Jo Ji Ji 



and the claim follows by Fubini's theorem. □ 
We note that the coefficient of t vanishes for k = 2. This jibes with the fact [3] 
that for any compact K in R m , m > 1, and all i > 

Q2,m(t) = 2<3l,m(t) — Ql,m(2t). 
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